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ABSTRACT 

We obtained over 90 hours of spectroscopic observations of four exoplane- 
tary systems with the Hobby- Eberly Telescope (HET). Observations were taken 
in transit and out of transit, and we analyzed the differenced spectra — i.e., the 
transmission spectra — to inspect it for absorption at the wavelengths of the neu- 
tral sodium (Na I) doublet at AA5889, 5895 and neutral potassium (K I) at A7698. 
We used the transmission spectrum at Ca I A6122 — which shows strong stellar 
absorption but is not an alkali metal resonance line that we expect to show 
significant absorption in these atmospheres — as a control line to examine our 
measurements for systematic errors. We use an empirical Monte Carlo method 
to quantity these systematic errors. In a reanalysis of the same dataset using a 
reduction and analysis pipeline that was derived independently, we confirm the 
previously seen Na I absorption in HD 189733b at a level of (— 5.26 ± 1.69) x 10~^ 
(the average value over a 12 A integration band to be consistent with previous 
authors). Additionally, we tentatively confirm the Na I absorption seen in HD 
209458b (independently by multiple authors) at a level of (-2.63 ± 0.81) x 10-^ 
though the interpretation is less clear. Furthermore, we find Na I absorption of 
(—3.16 ± 2.06) X 10""^ at < 3cr in HD 149026b; features apparent in the trans- 
mission spectrum are consistent with real absorption and indicate this may be 
a good target for future observations to confirm. No other results (Na I in HD 
147506b and Ca I and K I in all four targets) are significant to > 3cr, although 
we observe some features that we argue are primarily artifacts. 
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INTRODUCTION AND BACKGROUND 



Astronomers have now had evidence of planets orbiting other stars fo r roughly two 



decad es, including: the possibility of a planet around HD 114762 de scribed bvlLatham et al. 



(|l989[) ; evid ence of planets around pulsars due to timing variations ( Wolszczan fc Frail 



1992 



Rasidll994j); and the fa mous detection of a " hot Ju piter" around the solar-type star 51 Peg 



(IMayor fc Quelozl Il995l ) . The iLatham et al.l ( 1l989l ) and (IMayor fc Quelozl Il995l ) discoveries 
used the radial velocity (RV) method, observing the cyclic Doppler shifts in a star's spectrum 
to detect its orbit around a center of mass caused by one or more unseen companions. As of 
this writing, the RV method is still the dominant discovery method for confirmed detections 
of exoplanets. 

Over the last decade or so, astronomers have also observed the effects of exoplanets that 
pass in fro nt of their host stars and create measurable reductions in the observed light of 
those stars (jCharbonneau et al.ll2000l : iHenry et al.ll2000l : iMazeh et al.ll200Cll ). Transits are an 
extremely important phenomenon because they uniquely provide critical characteristics of 
planetary systems, such as the ratio of the planetary and stellar radii, and the ratio of the 
planetary orbit and the stellar radius. RV measurements, in contrast, cannot provide any 
information regarding the former, and can only be used to determine the latter if assumptions 
from stellar models are used. Planets that transit their central stars also present the only 
practical opportunity for astronomers to observe the properties of exoplanetary atmospheres 
in transmission. This is important because transmission spectroscopy of sufficient resolution 
is arguably the least ambiguous way to identify specific atomic or molecular species, thus 
providing important information about exoplanetary atmospheric composition. 

Many of the exoplanets that are known are "hot Jupiters," planets with masses ap- 
proximately that of Jupiter and with orbital semi-major axis a <C 1 AU. Such planets are 
observationally favored by searches using both radial velocity and transit methods. Con- 
sequently, these are also the planets most conducive to studies of their atmospheres using 
transmission spectroscopy. 

There are many different dimensions of the detection of an exoplanetary atmosphere. 
The central issue is the atomic or molecular species being detected. Detecting an atmosphere 
requires a differential measurement to distinguish the atmosphere from the signal of the disk, 
and thus all observations can be broadly classified as spectroscopy, ranging from narrow- 
band photometry to high-resolution spectroscopy. The wavelength, or more broadly the 
waveband (e.g. UV, optical, IR), dictates properties of the observation such as the telescope, 
instrument, and whether the observation is ground- or space-based. Detections may be 
made in absorption, using transits, or emission, using secondary eclipses and direct spectra. 
Finally, detections may be of a bound atmosphere at a low number of characteristic scale 
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heights or an unbound exosphere at many scale heights. All of these different dimensions are 
intertwined, but practically result in a very heterogeneous set of atmospheric observations 
that have been made at this early stage in the history of exoplanetary study. We discuss 
these results below. 



1.1. Overview of Existing Atmospheric Detections 



The literature contains a relatively small but rapidly increasing number of detections 
of a specific atomic or molecular species in an exoplanetary atmosphere. Due to obser- 
vational constraints , thes e detections have been primarily limited to only a few targets. 
Charbonneau et al.l ( 12002| ) was the first to detect an exoplanetary atmosphere. Using opti- 
cal spectroscopy from t he Space Telescope Imaging Spectrograph (STIS) onboard the Hubble 
Space Teles cope fffiSTl. |Charbonneau et al.l (120021) detected neutral sodium in HD 209458b. 



Soon after, IVidal-Madjar et al.l ( 12003 



20041 ) discovered the exosphere (the extended, un- 



bound portion of an atmosphere) of HD 20 9458b by detecting hydrogen , carbo n, and oxygen 
using the ultraviolet capabilities of STIS; iLecavelier Des Etangs et al.l ( 120101 ) made a sim- 
ilar detection of hydro gen in HD 189 7 33b u sing ACS spectroscopy. In contrast to these 
space-based detections, iRedfield et al.l (120081 ) made the first ground-based detection of an 



exoplanetary atmosphere us ing the Hobby - 



jber 



y Telescope (HET), detecting sodium in HD 
189733b. In the same year, ISnellen et al.l (l2008l) also made an im portant ground-based ob- 



servation, using Subaru to confirm the lCharbonneau et al.l (120021 ) detection of sodium in HD 
209458b. 

The IR has also been an extremely fruitful waveband for detecting exoplanetary atm- 
sopheres. For example, there have been detections of water vapor (H2O), carbon monox- 



i de (CO 



I Barman 



2010 



Tinetti et al. 



, carbon dioxide (CO^), and methane (CH4) in HP 209458b and HP 189733 b 
2OO7I: ISeaulieu et al.l 120081 boid ICrillmair et aPbood : ISwain et al.lbood . l2009bl 



20071 ). We note, how ever, that some o: 



are in significant dispute. For example, lEhrenreich et a. 



these detections that used Spitzer 
J20071) present an alternate inter- 



pretation of the Spitzer datas e t use d by iTinetti et al. (l2007l ) to claim a detection of H2O 
in HP 209458b. iGibson et al.l (120 111 ) also dispute molecular detections in HP 189733b and 
XO-lb that were made using HST/mCMOS. 

The above discussion clearly shows how the study of exoplanetary atmospheres has been 
focused on the two targets HP 189733b and HP 209458b. These systems are nearby and 
therefore the apparent brightness of their central stars enables transmission spectroscopy 
at greater signal-to-noise levels. More recently, however, there have been several additional 
detections of exoplanetary atmospheres. These detections include: observations of sodium in 
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WASP-17b flWood et al 



201 



2010f ) and XQ- 2b fISing et a 



of WASP- 12b flFossati et al. 



); the first detections of potassium, in HD 80606b (I Colon et al. 



20111) : and the presence of metals in the extended atmosphere 



20101 ). 



1.2. Alkali Metal Absorption in Exoplanetary Atmospheres 



With the detection of transiting exoplanets, e arly models ( ISeager fc Sasselovl 12000 



Barman et al.ll200ll : lBrowrul200ll : iHubbard et al.ll200ll ) sought to predict what features might 
dominate the atmospheres of hot Jupiter-type exoplanets. In the optical, the two dominant 
features should be t h e reso nance absorption lines of the alkali metals sodium and potassium. 
Charbonneau et al.l ( 12002| ) sought to test this prediction by observing the Na D doublet 
at AA5889, 5895. They used STIS at a spectroscopic resolution of R = X/6X = 5540, and 
created a transmission spectrum of the atmosphere of HD 20945 8b by subtracting th e coad- 
ded in-transit spectra from the coadded out-of-transit spectra. ICharbonneau et al.l (120021 ) 
found an absorption level of (2.32 ± 0.57) x 10~^ integrated over a 12 A band relative to 
the adjacent spectral regions, with dec reasing levels of abs orption detected as the band- 
width was increased. From the ground, [Snellen et al.l ( 12008! ) confirmed this result, finding 
a (5.6 ± 0.7) X 10~^ absorption level summed over two 3 A bands centered on the two 
lines of the Na D doublet , which ISnellen et al.l (120081 ) determined was consiste nt with the 
Charbonneau et al. ( 20021) result when differences in bandpasses are considered. Sing et al. 



(l2008al ) reanalyzed the ICharbonneau et al.l (120021 ) data and calculated the absorption level 
over several different bandpasses. They concluded that bot h narrow and broad absorption 
was present, at a level cons istent with early cloudless models (ISeager fc Sasselovll2000l : iBrown 



2001 



Hubbard et al.ll200l[ ). 



Redfield et al.l (120081 ) detected sodium in the atmosphere of HD 189733b, at a level 
of (6.72 ± 2.07) X 10-"^. This measurement is an integration of the transmission spectrum 
5*^ = {Fin/ Font) — 1, where Fin and Font are, res pectively, the r iorma lized in- and out- 



of-transits fluxes. The absorption level detected by iRedfield et al. J2008[l was recent l y con- 
firmed by i/S'T measurements (D. Sing, 2010, private commu nication). iRedfield et al.l (l2008r) 
integrated St over a 12 A bandpass for consistency with the ICharbonneau et al.l ( 20021 ) ob- 
servations of HD 209458b. We do not necessarily expect these atmospheres to be similar, 
as infrar ed observations indic ate that the atmosphere of HD 209458b has a temperature in 



versio n ( iKnutson et al.ll2008l ) but t he atmosphe r e of H D 189733b does not (iGrillmair et al. 



20071 ). It is important to note that iPont et al.l (120081 ) found no distinct features using the 
Advanced Camera for Surveys (ACS) onboard HST, implying a broad haze that does not 
allow for an absorption signature. However, the low spectral resolution {R ~100, compared 
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to R ~60,000 in lRedfield et al.l 120081 ) might mask the absorption 



A third detection of sodium was made recently by IWood et al.l (1201 ll ). who found ab- 
sorption in WASP- 17b at a level of (5.5±1. 3) x 10~^ using a 2 x 1.5 A bin. This is significantly 
less than predicted by taking the model of iBrowru (120011 ) applied to HD 209458b and scaling 
it by the difference in scale height of the two planets. Also notable is that the absorption is 
significa nt only using narro w bandpasses; as the bandpass is increased, the absorption dis- 



appears. 



Wood et al. 



(120111 ) discuss several possibilities to explain the discrepancy between 
models and their observations, including photoionization on the day side of the planet, the 
depletion of heavy elements on the cold night side, and clouds. In particular, clouds could 
explain the nature of the absorption as a function of bandpass width, as broad absorption 
should be formed at high pressures and low altitudes and narrower absorption at lower pres- 
sures and higher altitudes. Thus, a layer of clouds or other haze at the proper height would 
suppress the broad absorption but not the narrower absorption. 



Recentl y, the first detect ions of potassium were made for HD 80606b ( IColon et al.ll2010l ) 
and X0-2b (jSing et al.ll201ll ). Both detections where made usi ng tunable- fi lter p hotometry 
with the OSIRIS instrument on the Gran Telescopio Canarias. ISing et al.l (120111 ) find clear 
absorption in the K I core at 7664.9 A and in the blue K I wing at 7582 A. They find 
the a pparent plarietary radius to be consistent with cloudless models for X0-2b. In con- 



trast, IColon et al.l (120 lOf ) see broad absorption but no additional absorption in the K I core, 
implying that low-pressure, high-altitude potassium has been ionized or has condensed. 

Taken as a whole — and given that broader searches in the optical have not revealed any 
other significant features — these results confirm that neutral sodium and potassium provide 
the dominant line opacity sources in the optical, in accordance with models. They also 
highlight the issues that still exist. In a given planet, what is the effect of photoionization 
and condensation on abundances? How common are clouds and haze? New observations 
of additional targets, preferably at high resolution, will be required to fully address these 
questions. 



1.3. Future Prospects and This Study 

The trends in the detection of exoplanetary atmospheres discussed in §l.ll and §1.2l show 
that we are entering an exciting time for this subfield. A variety of techniques have now been 
demonstrated to be feasible (ground- vs. spaced-based, photometric vs. spectroscopic, etc.). 
Having an array of techniques and an ever-increasing target list means we are on the cusp of 
performing significant comparative atmospheric exoplanetology for the first time. It would 
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be particularly useful to have a homogeneous survey of many different targets in order to 
compare absorption in exoplanetary atmospheres, and determine if variations are correlated 
with stellar or planetary characteristics. 

The study described in this paper aims to be a first step toward such a survey. We 
have made uniform observations of the spectra of four exoplanetary systems. The goal is to 
examine, with high-resolution spectroscopy, the transmission spectra of several exoplanetary 
atmospheres and search for absorption from the Na I doublet at AA5889, 5995, along with the 
K I absorption line at A7698|l| We describe the observations and basic data reduction in §2J 
In ^we describe our procedures for analyzing the reduced spectra, along with a discussion 
of vetting our results, including error analysis and examining "control lines" for comparison. 
In §1] we describe our basic results. We discuss these results in §S] and summarize in §Si 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. Description of Observations 

Observations were obtained with the HET between August 2006 and November 2008. 
The driving factor behind our use of the HET is its large primary mirror with an effective 
aperture of 9.2m. This is large enough to provide the high signal-to-noise that we require 
in these observations. Spectra were observed with the High- Resolution Spectrograph (HRS) 
at a nominal resolution of i? ~ 60, 000, though we measure wavelength-dependent and time- 
variable resolutions at our specific wavelengths of interest (see §2.2p . Individual exposures 
were limited to 10 minutes in length in order to avoid saturation, though some are shorter. 
Approximately 20% of the exposures for each target were taken while its respective planet 
was in transit. Specific details are provided in Tables [1] (observational statistics) and [2] 
(system and orbital parameters). 

We have an interest in spectral features that are near telluric lines and we also have 
very high signal-to-noise requirements. Therefore, part of our observing strategy was to 
observe "custom" telluric standard stars to use in removal of Earth's atmospheric lines. 
These stars were chosen to be near our exoplanetary targets on the sky and the observations 
were taken either immediately before or after our primary science observations. The spatial 
and temporal proximity reduces differences in air mass and water vapor content between the 
primary observations and the telluric observations. 



^Note that the K I Hne is also part of a widely-spaced doublet, but the other component at A7664 is close 
to much stronger telluric absorption than the 7698 A line. 
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2.2. Data Reduction 

The raw echelle spectra from the HET were reduced using standard IRAF processes to 
remove the bias level of the CCD, flatten the field, remove scattered light, extract the spectra 
from the apertures, and determine the wavelength calibration from the ThAr lamp exposure. 
The primary science targets, telluric observations, and ThAr exposures were all extracted to 
one-dimensional spectra. We show examples of reduced but pre-processed primary science 
and telluric standard spectra in Fig. [H The examples shown are for HD 189733 and its 
telluric standard HD 196867; the entire echelle orders containing Na I, Ca I, and K I are 
shown. 

Before subtracting the telluric standard observations from the primary science observa- 
tions, we must "clean" them of weak stellar and interstellar lines. We fit simple models to the 
extracted telluric star observations, simultaneously modeling the star and the tell uric absorp- 



tion lines using IDL programs. These programs used the MPFIT architect ure (iMarkwardt 



200 



^ for minimization. The star is modeled with a b ase stellar model ( Koesterke et al 



20081 : iFitzpatrick fc Massall2005l : ICastelli fc Kuruczll2006l ) with free parameters of a velocity 



shift and rotational parameter vsini. Absorption from the Earth's atmosphere is modeled 
with a base telluric model with free parameters of a velocity shift and optical depth. We ex- 
amine the residuals of these fits for additional absorption. Among our three telluric stars, we 
find various combinations of interstellar absorption and additional stellar absorption (likely 
from a companion to the telluric standard star), with wavelengths corresponding to the Na 
D doublet and K I A7698. These lines are fit for column density and 6- value (i.e., Doppler 
parameter) in order to characterize their profile imprint. 

For a given telluric star, all observations are fit individually at first, with all free pa- 
rameters described in the previous paragraph allowed to vary. Then the median fits of v sin i 
for the star and the median column density and 6-value combination of the additional inter- 
stellar and stellar companion components are selected. These values are then fixed and the 
observations are re-fit for velocity shifts and telluric optical depth. This new, standardized 
best fit of the stellar model and interstellar /stellar companion absorption are then divided 
from the telluric observation. The telluric absorption is modeled in order to aid the fit but is 
not applied to this division, resulting in a telluric observation "cleaned" of stellar and ISM 
absorption imprints. 

These "cleaned" telluric observations are then removed from the primary science ob- 
servations using standard IRAF procedures. In cases where our custom telluric star was 
not observed in a given night, we compare the primary observation to all available telluric 
observations of the same custom telluric star over a region dominated by telluric lines and 
select the match with the lowest (after a fit for velocity shift and relative flux scale). 
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Removing the telluric lines removes the majority of the imprint of the blaze function 
on the primary science observations. We then use a low-order polynomial to complete the 
normalization and fix the scale to unity. Flux errors are calculated based on photon and 
read noise, and scaled by the normalization factors that are introduced from the telluric 
subtraction and additional normalization. 

We coadd all observations into two sets, in-transit and out-of-transit. Based on the 
ephemeris we use (see Table H]), some observations cross through first or fourth contact 
during the exposure. We classify these observations as either in- or out-of-transit based on 
the midpoint of each observation. Observations with a midpoint during ingress or egress are 
classified as "in-transit," although we make note of how many observations this applies to 
in Table [U We also test whether or not using only observations with a midpoint during full 
transit (i.e., between second and third contact) makes a difference in the final results (see 

During our coadd, we apply two higher-order corrections. First, we observe that the 
resolution of the HET at our settings varies with time and with wavelength. This variation 
is most severe at short wavelengths of a given order, with the greatest variation before 
December 2007. We empirically measure the resolution by fitting the width of unsaturated 
lines in the ThAr lamp exposure. We measure lines over the entire wavelength region of each 
order and then fit the resolution as a function of wavelength in each order for each night with 
a low-order polynomial. At each resonance line of interest, we find the minimum resolution 
over all observations and degrade the remaining observations to a chosen resolution slightly 
less than this using a Gaussian kernel (e.g., we find the resolution of all observations to be 
at least slightly greater than R = 65, 000 near Na D, and thus we degrade all observations 
to R = 65,000). 

The second correction that we apply is a wavelength correction between the various 
observations. This corrects for both any overall shift in the observations and higher-order 
variations in the wavelength solution that can be observed by inspecting the stellar lines. 
We perform a cross-correlation on strong stellar lines and fit a spline through the cross- 
correlation values. This correction is applied to the wavelength arrays and the fiux and error 
arrays are interpolated onto the new wavelength grids. The final coadding is weighted by 
the errors at each point in the spectrum, with the errors based on photon statistics and read 
noise, as described above. 
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3. ANALYSIS 

3.1. Transmission Spectra 

As described in the previous section, we make two coadds — one of all in-transit spectra, 
and one of all out-of-transit spectra. We perform another wavelength correction as described 
in §2.2l in order to ensure that the two wavelength arrays are on identical grids and correct for 
any wavelength solution differences that still exist between the two sets of coadded spectra. 
We then divide the in-transit spectrum by the out-of-transit spectrum and subtract one 
(to shift from a spectrum normalized to unity to one normalized to zero) to calculate our 
"differenced" or "transmission" spectrum. 

We averaged the transmission spectra over several different fixed wavelength regions 
or "wavebands" — a central waveband centered around the line we want to analyze and two 
additional wavebands of equal size, immediately to the blue and red of the central waveband. 
The absorption lines and the wavebands widths are given in Table [31 The center of the 
waveband is determined from the minimum point of the absorption line, which varies slightly 
between our targets and when examining subsets of the data (see §3.2p : in the case of the 
Na I doublet, the midpoint of the waveband is the midpoint between the minima of the two 
doublet lines. The total absorption level is then the central waveband minus the average of 
the blue and red wavebands. We carry through, using standard error propagation, the formal 
statistical errors on each point through the differencing and averaging steps to calculate a 
statistical error of the total absorption level. 



3.2. Estimating Errors with Monte Carlo Methods 



Our final results are sensitive to any systematic errors that may be present in the reduc- 
tion process — extraction of the spectra, cleaning of the telluric standard spectra, removing 
the telluric imprint from the primary science spe ctra, coadding the spe ctra, and differencing 
the coadded in-transit vs. out-of-transit spectra. iRedfield et al.l (120081 ) suspected that these 
possible systematic e rrors are larger than the formally propagated statistical errors. Thus, 
Redfield et al.l (120081 ) devised an empirical check on systematic errors. They referred to their 
method as an "Empirical Monte Carlo" method (hereafter EMC), and we repeat this method 
of analysis here. 

We use three versions of the EMC analysis. In each version, we perform the analysis 
described in §3. II but instead of the total coadded in-transit spectrum, a subset of observa- 
tions is taken, coadded, and a differenced spectrum is calculated relative to a single "master" 
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spectrum. We perform thousands of iterations of this process for each EMC method in order 
to examine the centroid and width of the distribution and see how this compares with the 
master transmission spectrum of each target and absorption hne combination. 

The first EMC method we refer to as the "in-out" method, where we randomly select 
a subset of in-transit spectra, coadd the selections using the same methods as above, and 
then calculate a transmission spectrum relative to the total coadded out-of-transit spectrum. 
The number of selected spectra in the in-transit subset is randomly varied, from five up to 
one-half of the total number of in-transit spectra available, e.g., HD 189733 has 35 in-transit 
and 172 out-of-transit observations (observations which are not discarded; see Tabled]), so 
between 5 and 18 in-transit spectra are chosen for the subset. We expect the in-out method 
to result in a distribution with a centroid consistent with the absorption level calculated 
including all observations. 

The second method is the "out-out" method, where the calculated transmission spec- 
trum is not a true transmission spectrum, but a comparison of a randomly selected subset of 
out-of-transit observations, compared to the complete coadded out-of-transit spectrum. The 
number of observations in the subset is the same as the number of total in-transit observa- 
tions (e.g. 35 for HD 189733). In principle this method should show a normal distribution 
centered at (no absorption) with a width that indicates a robust measurement of systematic 
errors in our data processing and the system itself (e.g. stellar variability). 

The final method is the "in-in" method. As with the out-out method, a random subset 
of in-transit observations are compared to the complete in-transit spectrum. The number 
in the subset is chosen such that the number of subset observations to total number in- 
transit observations is the same as the ratio of the number in-transit observations to the 
number of out-of-transit observations (e.g., 7/35 ~ 35/172 for HD 189733). As with the 
out-out method, the in-in method should be centered at no absorption. The width of the 
distribution, however, is less meaningful because of the lower signal-to-noise level resulting 
from the smaller number of observations. 

This trifecta of methods is differe nt from the more common "bootstrap" or "quick-and- 



dirty" Monte Carlo (iPress et al.lll992l ). In the bootstrap Monte Carlo, a dataset of an equal 
number of points as the original set is chosen, with replacement possible. "Replacement" 
means that a randomly chosen point (or in our case, observation) is "put back (replaced) 
into the bin" from which the points are drawn, and therefore duplication is possible, and 
in a typical iteration 1/e ~ 37% points in the dataset will be duplicated. The methods 
described above use randomly chosen subsets without repla cement, and therefore bear some 



similarities to a "jackknife" method ( iWall fc Jenkins! |2003| ). There are two primary reasons 



that we use our "out-out" EMC method to determine the error. The first is that by using 
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the "out-out" method we are exphcitly removing the effects of planetary transit when we 
attempt to quantify the astrophysical noise (in addition to the statistical errors and reduction 
systematics). The second reason is that by using the "out-out" method, we are able to draw 
upon a larger set of points for the Monte Carlo analysis. However, we also test the more 
traditional bootstrap Monte Carlo, in a way that is most analogous to our "in-out" method. 
Using HD 189733b as an example, the boostrap method would entail a random sampling of 
35 "in-transit" observations, with replacement (i.e., duplication) allowed, and a sample of 
172 "out-of-transit" observations, again with replacement. We discuss these tests along with 
the rest of the results in §H 



3.3. Confirming Detections with a Control Line 



Redfield et al.l ( l2008l ) also tested the potential objection that there may be systematic 
errors simply from the process of coadding and differencing spectra surrounding strong stellar 
absorption lines. In order to test this, they calculated a transmission spectrum at Ca I A6122, 
which is a reasonably strong stellar line in our target systems but not a resonance line that 
is expected to manifest in hot Jupiter-type atmospheres according to models (it is expected 
to be condensed out, e.g., into CaTiOa). We also calculate a transmission spectrum of this 
line for all four targets. 

The wavelength region near this line is utterly devoid of telluric features, so our data 
analysis process is modified slightly in that we do not perform telluric subtraction. This 
is advantageous in that suboptimal telluric subtraction will not influence the results. It is 
disadvantageous in that it means the control line analysis is different from the Na and K 
line analysis in this particular step (all other steps are identical). However, it provides a 
consistency check on the quality of our telluric subtraction. 

In addition to calculating a total transmission spectrum, we also run the EMC analysis 
for this line. We expect a null master measurement (no absorption) and an EMC distribution 
centered around zero as well to confirm that our measurements are not simply a typical 
systematic byproduct of differencing around strong stellar lines. 



4. RESULTS 

We have four targets that we examine for evidence of Na I (AA5889,5895), Ca I (A6122), 
and K I (A7698). We expect Na I and K I to exist in atomic form in hot Jupiter atmospheres 
and for these ground-state resonance lines to exhibit absorption. However, as described in 
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the previous section ( §3.3p . this is not true of Ca I, which functions as a control hne. The 
transmission spectra for each hne and target are shown in Fig. [SHU These are shown in order 
of increasing wavelength: Na I, Ca I, and K I, respectively. Our measurements are shown in 
Table |2 

The corresponding EMC error analyses are shown, in the same order, in Fig. |5HZ1 In 
the discussion of the subsequent sections, we will use the width of the "out-out" EMC 
distribution as the la error. We first discuss our results for Ca I, the control line. 

4.1. Ca I Control Line 

As discussed in §3.3[ we do not expect absorption due to Ca I, as calcium should be 
condensed out at the atmospheric temperatures of our targets. In three of the four targets, 
the master measurement of flux at the control line is consistent with zero to within the 
la error, as we expect. In HD 149026b, the master measurement deviates from zero, at 
a significance level of lA8a. The deviation is positive, indicating "emission" as opposed 
to absorption, but such a result, if robust, is not physically meaningful for transits. Minor 
features can be seen in the top right panel of Fig. [3l consistent with the calculated significance 
level. 

Using a standard bootstrap Monte Carlo method provides similar results — 3 of the 4 
targets are within la of zero, and HD 149026b is within 2a. Notably, the error (i.e., the 
distribution width) produced by the bootstrap method is slightly larger than the "out-out" 
EMC width, but we have explained in §3.21 our reasons for adopting the latter as our error. 
In addition, if we perform a master transmission measurement using only observations with 
a midpoint between second and third contact, the results are again essentially the same. The 
results for all four targets taken as a whole confirm that statistically significant absorption 
and emission profiles are not inevitable, or even probable, simply as a byproduct of our 
reduction and analysis processes. This provides confidence in the results that we do detect 
in other lines. 



4.2. Nal 

We detect absorption at Na I in HD 189733b at a significance of 3. Oct, with an absorption 
level of (-5.26±1.69) x 10"^ This is consistent to within 0.55a with the (-6.72±2.07) x 10"^ 
absorption found by R2008 and confirmed by recent HST measurements (D. Sing, 2010, 
private communication; Huitson, et al., submitted). The profiles can be seen in the bottom 
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left panel of Fig. [2] and are comparable to the profiles seen in R2008. We note that our 
original data is the same as that used by R2008. However, our data reduction and analysis 
were not identical, using a quasi-independently derived set of IRAF reduction scripts and 
IDL analysis programs in order to apply these programs to the broader dataset of additional 
targets and K I analysis. 



^HD 20 9458b is also known t o ^how Na I absorption ( ICharbonneau et al.ll2002l : [Snellen et al, 

2008 ). The Charbonneau et"aD (2002) detection was (—2.32 ±0.57) x 10""^ in an absorption 



band identical to ours, while [Snellen et al.l (120081 ) used a variety of bandpasses that are not 
directly comparable to our results. We fin d abso rption at a level of (—2.63 ± 0.81) x 10""^, 
nearly identical to the ICharbonneau et al.l ( 120021 ) measurement. However, the transmission 
spectrum (bottom right panel of Fig. [2]) is less clear, with only the narrowest features seen 
at the Na I lines, along with several other narrow features in the red comparison band. In 
addition, of the 12 EMC analyses (4 targets and 3 absorption lines), this is the only case 
where the centroid of the "in-out" distribution is discrepant with the master measurement, 
as the "in-out" centroid is consistent with zero ((—0.57 ± 0.81) x 10~^). A possible ex- 
planation is the narrow features in the red comparison band, as our master measurement 
is based on the integration of the central bands defined in Table [3] minus the average of 
equivalent-width bands to immediately to the blue and red. If we modify the calculation so 
that we only subtract the blue comparison band for both the master measurement and the 
EMC results, then the centroid of the "in-out" distribution changes to —1.44 x 10"^ and 
the master measurement to —3.30 x 10^^. These results are 0.9a and l.Oo" discrepant from 



the ICharbonneau et al.l (120021 ) results, respectively, but still discrepant from each other. It 
is possible that the narrow features in the red comparison band are artifacts of the telluric 
subtraction process, but this is not clear. The three strongest of these features roughly match 
known telluric lines, but there are inconsistencies in the velocity offsets. Furthermore, we 
do not see features in the transmission spectrum at the wavelength of other telluric lines 
of nearly comparable strength. We have also examined the individual spectra and are un- 
able to determine if and in what spectrum the possible subpar telluric subtraction occurs. 
Therefore, it is unclear if poor telluric subtraction could be responsible for these features. 
Currently we do not have a clear explanation for the cause of these features or the internal 
inconsistency between the master measurement and the EMC centroid. 

Our measurement for HD 147506b indicates a 1.6cr surplus. However, no feature is 
apparent in the transmission spectrum (top left panel of Fig. [2]). Conversely, HD 149026b 
seems to show absorption at a comparable level of 1.5a. Though the significance levels 
are similar, the latter has a greater magnitude and a larger error. Additionally, absorption 
features are seen at both Na I lines (top right panel of Fig. [2]). We hypothesize that this 
absorption may be real, though additional observations would be required to confirm this. 
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Similar to the results for Ca I, the bootstrap Monte Carlo results for Na I are consistent 
with our master and EMC results. In the case of HD 209458b, the bootstrap distribution is 
like our "in-out" EMC distribution in that it is consistent with zero and inconsistent with 
the master measurement. The bootstrap errors are also slightly larger (< 50% on average) 
than the "out-out" EMC results. 

If ingress and egress observations are discarded, the results are still consistent with our 
original master measurements to within <lcr. We do note that the measured HD 189733b 
Na I absorption increases to —6.50 x 10~^. This is notable in part because of our four 
targets HD 189733b has the greatest percentage of ingress/egress observations that we must 
discard (due to the planet's larger Rp/R* and shorter orbital period, the ratio of time spent 
in ingress/egress to full transit is larger). Thus, if signal dilution from the ingress/egress 
observations is a problem for our original master measurement, we would expect the effect 
to be strongest in this case. The Na I absorption levels for HD 149026b and HD 209458b, 
however, decrease to —2.04 x 10"'^ and —1.87 x 10~^, respectively. 



4.3. KI 



(Colon et al. 


2010; 


Sine et al. 


2011) 



absorption, though our measurement for HD 209458b is (—0.89 ±0.67) x 10 ^, which is 1.3cr 
significant. HD 189733b shows no feature and is consistent with zero to within the la errors. 

The other two targets, however, do show features worth mentioning. HD 149026b show 
a feature in the blue comparison band that is visually significant relative to the noise. This 
feature is not centered on any feature in the stellar spectrum, and is too far from K I in 
velocity space to be blueshifted K I absorption. If we neglect the blue band and compare 
the central measurement to only the red band, then the master measurement of the K I is 
0.02 X 10~^, miniscule compared to the 1.09 x 10~^ error. The feature itself, integrated over 
the entire blue band, is also less significant than the la error. 

HD 147506b is the most difficult to understand case. There is an emission-like feature 
centered on K I in the transmission spectrum (upper left panel. Fig. Hj). According to our 
EMC calculations, this feature is 2.1a significant. There is no easily understood physical 
reason why a transmission spectrum should show apparent emission. However, we note 
that of our four targets, HD 147506b has relatively low overall S/N in the raw observations 
(comparable to HD 149026b but much lower than HD 189733b and HD 209458b). Its stellar 
features are also the broadest and shallowest of our entire sample (due to its v sini of 20.8 ± 
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0.3 km s~^, > 3x the other three stars in our sample; |Pal et aLll2010l : iButler et al.ll2006l ). 
which comphcates the hne cross-correlation required for the wavelength corrections we use to 
match up the master in- and out-of-transit spectra before performing the differencing step. 
Furthermore, the echelle order containing K I for HD 147506b was by far the most uncertain 
dataset qualitatively, as we discarded multiple observations due to reduction artifacts (spikes 
and dips in the nearby spectrum) of unknown origin. Though these were discarded, it is 
unclear if any observations remain that unduly influence the measurement here. For these 
reasons, and the fact that feature is only 2.1a significant, we argue that this feature does 
not correspond to any real physical effect, and is most likely isolated to the analysis of this 
particular absorption line in this single target. 

Once again, considering the bootstrap Monte Carlo or only the transits only between 
second and third contact do not appreciably change any of the results. We note, however, that 
the much larger error from the bootstrap method would make the HD 147506b "emission" 
only 1.3(7. 



5. DISCUSSION 

We synthesize our results for all three elements in Fig. |Hi In this plot, it is easy to see 
the consistency of Ca I with zero, which we expect as our control line. The other results 
discussed in §4.2H4.3l can also be seen, including the > 3a detections of Na I in HD 189733b 
and HD 209458b and the very tentative {1.5a) indication of Na I absorption in HD 149026b. 
Fig. [9] shows the same data from a different perspective; it is similar to Fig. [HI except that 
plotted on the y-axis is the apparent radius of the planet, averaged over the bandpasses 
defined in Table [31 Fig. M also includes integrated model values, which we discuss later in 
§531 



5.1. Dependence on Scale Height 

Scale height is used in Fig. [HI and M because it is a straightforward measurement that 
should correlate with the amount of observed absorption, all other factors being equal. It 
is obviously difficult to pull out any clear trends from these two plots given that only two 
points are statistically different than zero (> 3a). However, we can say that our detections 
(Na I in HD 189733b and HD 209458b) and our nominal detection (Na I in HD 149026b) are 
consistent with what we expect from consideration of atmospheric scale height. We observe 
nothing for HD 147506b, which is the most massive of our four planets (but similar in size to 
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HD 189733b and HD 209458b), thus giving it the largest surface gravity and smallest scale 
height, while the other three planets, with much larger scale heights, show at least the hint 
of absorption. 

There are other properties beyond scale height that determine the detectability of at- 
mospheric absorption. One important property is the metal abundance. HD 149026 has 
the highest metallicity, [Fe/H] = 0.36 ± 0.08, of all the central stars of our targets. HD 
149026 b can likewise be inferred to be particularly metal rich based on bulk density argu- 
ments ( ISato et al.l 120051 ). While there is no guarantee that these quantities would translate 
into an enhanced metal abundances in the atmosphere, bu lk density and atmospheric abun- 
dances are correlated in the Solar System ( iLoddersI 120031 ). It is therefore possible that HD 
149026b might have an enhanced atmospheric Na I abundance, which could be responsible 
for the large (albeit statistically tentative) R^/Rp that we derive for the 12 A wide band 
around Na I in this target. 



5.2. Temperature Inversions 



Another important consideration is the issue of atmospheric temperature inversions. 
Spitzer observations of many extrasolar planets have revealed two classes of planets: those 
whose IR spectra are fit well by models with H2O and CO in absorption, and planets 
whose spectra are better fit by models with th ose same molecules in emission, indicat- 
i ng a temperature inversio n (IHubeny et al.l 120031 ). HD 189733b is i n the former categor y 
jCharbonneau et al.l mKM ) and HD 209458b is in the latter class jKnutson et al.l I2OO8I ). 
There has been comparatively little IR work on HD 14750 6b and HD 149026b, with no 
Spitzer observations of the fo r mer (ISeaeer fc Deming I2OIOI ) and on ly 8 nm work on the 



latter ( Harrington et al. 2007; Knutson et al. 20091 ). Consequently, Knutson et al. ( 2010 ) 



classifies these planets as unknown with respect to temperature inversions. 

A temperature inversion is relevant to Na I and K I absorption because it requires an 
additional absorber ( e.g. VO/TiO or clou ds ) that may obs c ure th e abs orption of other ma- 
jor opacity sources. iFortney et al.l (120031 ) . iBurrows et al.l (120071 ) . and iDesert et al.l (120081 ) 
have invoked such an absorbe r as a possible reason w hy the absorption in Na I in HD 
209458b is less than expected (jCharbonneau et al.ll2002l ) and less than that of HD 189733b 
( iRedfield et al.ll2008l ). which does not appear to have an inversion. If this is the case, then 
our possible detection of Na I in HD 149026b — with a scale height greater than HD 189733b 
but much less than HD 209458b — might nominally predict that it does not have an inversion. 



Knutson et al.l (120101 ) has proposed that increased stellar UV flux may destroy whatever un- 
known material is responsible for the inversion. Notably, the central star HD 149026 has 
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a smaller UV flux more consistent with the inversion cases shown in iKnutson et al.l ( 120101 ) . 
However, ISeager &: Deming ( 120101 ) note that there are counterexamples of high stellar activ- 
ity and assumed inversions, most notably XO-lb. They also note that many cases of inferred 
temperature inversions have adequate alternative explanations to within the uncertainties, 
e.g., an increased CH4 abundance relative to CO, and that three-dimensional circulation 
models are inadequate to explain all of the Spitzer/IRAC data, even with inversions in- 
cluded (see lSeager fc Demingll2010l . and references therein). We also note that other factors 
likely play a role in depleting th e observ ed at mospheric abundanc e s, inc luding ionization and 
condensation. Both ISing et al.l ( l2008bl ) and IVidal-Madjar et al.l ( 1201 ll ) conclude that there 
must be either an ionization or condensation effect in HD 209458b; in particular, the latter 
invokes both explanations at different layers in the atmosphere. 



5.3. Exospheric Absorption 



We have assumed that any detection of Na I or K I would come from the bound atmo- 
sphere of the planet. Measuring wavelength-dependent apparent planetary radius (i.e., Ra) 
larger than the Roche radius would indicate that some absorptio n comes from the unbound 



exosp here, e.g. the hydrogen detection in HD 209458b made by IVidal-Madjar et al.l ( 12003 



20041 ) through Lya measurements. However, the converse is not true; if Ra lies inside the 
Roche radius, the absorption may be either bound or unbound depending on the details of the 
absorption (e.g., abundance, saturation, etc.). Based on our transmission spectra in Fig. |2H1] 
and averaged results in Fig. IH we do not flnd evidence that any of the absorption is from the 
exosphere of the planet. However, given that both HD 189733b and HD 209458b appear to 
ha ve some escaping atmosphere, and are possibly in a state of "blow- off' (detected in Lya 
by Ihecavelier Des Etangs et al.l boiol : Ividal-Madjar et all I2OO3I . |2004| . respectively, for HD 
189733b and HD 209458b), it is possible that some Na I absorption occurs in the exosphere 
due to the effects of hydrodynamic drag. 



5.4. Velocity Shifts of Atmospheric Lines 



Redfleld et al.l (120081 ) detected a substantial apparent velocity shift of —38 km s ^ in 



the Na I lines of HD 189733b, relative to the Na I lines of the star. The authors suggested 
that the shift, if astrophysical in origin, may have been a combination of planetary motion 
and winds from the hot dayside fl owing to the nightsi de. More recently, such shifts have 
been measured in HD 209458b by ISnellen et al.l (|2010[ ) , although at a much smaller level 
of 2 km s~^. In this analysis, we measure apparent shifts in the lines by fltting Gaussian 
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profiles to the features we noted near Na I in HD 149026b, HD 189733b, and HD 209458b 
(altliougli only tlie latter two are detected at significant levels; see §4.2p . We find shifts 
of approximately —4, —37, and —12 km s~^, respecti vely. In the case of HD 189733b, our 
result is identical to the result of iRedfield et al.l ( 120081 ) to within the precision of the spectral 
resolution, which is ~ 5 km s"^ for the HRS at its R - 60, 000 setting. 

There is a qualitative difference between HD 189733b and the other two targets, in 
that the feature we associate with Na I absor ption is highly asymrn etric. This asymmetry is 
seen both in the reduction of this paper and iRedfield et al.l (120081 ). Before performing the 
step of differencing the in- and out-of-transit spectra, we perform wavelength corrections by 
performing cross- correlations on the lines between the two spectra to create a wavelength 
correction array. If this cross-correlation is imperfect, then there will be an asymmetric 
feature created, with a dip to one side of the stellar line center and a spike to the other side, 
which is what is seen in our data (with a dip to the blue and a spike to the red, though the 
integration yields net absorption). Whether or not this is the source of the artifact we see, 
we note that fitting a Gaussian to an apparent absorption feature when this type of artifact 
is present will artificially force the fit of line center to the side with the dip. In HD 189733b, 
this is to the blue, whereas HD 149026b and HD 209458b do not show the same type of 
asymmetric artifact. If instead of fitting Gaussian profiles we assume that line center of the 
atmospheric feature is the local minimum point in the transmission spectrum, then we find 
shifts of approximately 0, 1, and -11 km s"^ to HD 149026b, HD 189733b, and HD 209458b, 
respectively. The results for HD 149026b and HD 209458b are consistent with the Gaussian 
results to within the ~ 5 km s~^ resolution element. We conclude that it is the nature of 



the HD 189733b artifact, seen both in this present work and in IRedfield et al.l (120081 ). that 
causes the large apparent shift when measured through Gaussian line fits. This suggests that 
the present data may be inadequate to robustly make precise measurements of atmospheric 
velocity shifts, especially when an asymmetric artifact of the differencing step is seen. 



5.5. Comparison with Models 



In Fig. [9], there are additional points that show the integrated points of planetary atmo- 
sphere models. One model includes absorption by TiO and VO (which we call the "primary" 
model), while the other does not (the "no TiO/VO" model). The primary models, by their 
inclusion of these molecules, are imposing an absorber in order to creat e an atmospheric tern - 



perature inversion, which has been observed for several planets (e.g. iKnutson et al.l l2008l ) . 
We note, however, that there is ongoing debate about whet her or not TiO/VO can create 
a temperature inversion in a planetary amotsphere (see, e.g.. ISpiegel et al.ll2009l . who argue 
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the negative case). 

For HD 147506b, it can be seen that the model points are all very close to 1 in R;j/Rp, 
i.e., no additional absorption from the atmosphere. We noted in §l]that the integration of 
the transmission spectrum for this star gave an excess at Na I and K I, though less than 3a 
in each case. To the extent that these results are consistent with zero (in terms of relative 
absorption; 1 in R^/Rp), they are also consistent with both the primary and no TiO/VO 
models. We note that HD 147506b was included in this study because of its bright central 
star rather than planetary properties that suggested potential absorption]^ The fact that we 
see no obvious absorption is consistent with its much smaller scale height compared to our 
other targets. 

HD 149026b shows the most significant spread in R^/Rp, because its planetary radius 
is the smallest of our sample. The primary model values are much closer to zero, but 
consistent with our measurements within 1.5a for Na I, l.Qa for Ca I, and 0.65a for K I. 
The model without TiO and VO provides a slightly better fit to both Na I (within l.lcr 
vs. 1.5a) and K I (within O.Ola vs. 0.650"). We have suggested a possible detection of Na I, 
significant at the 1.5cr level but with a profile that suggests real absorption, and suggest 
that further observations should be performed to confirm or rule out this measurement. Our 
current measurement implies a case closer to the model without TiO/VO, though it does 
not strongly differentiate between the two models. 

The primary and no TiO/VO models are virtually identical for HD 189733b, because the 
temperature of this planet precludes the existence of these molecules. The models predict 
less Na I and more K I absorption than we see — 1.3cr and 2.7a differences. Despite the 
overprediction of K I relative to our measurement, the magnitude of our prediction implies 
it is unlikely we would be able to detect any absorption at significant levels. 

HD 209458b shows a good match at Na I (within 0.68cr) and Ca I (within O.Ala) and 
a fair match at K I (within 0.94cr) for the primary model. The model without TiO/VO 
does much worse, as it is 3.6ct and 4.5(j discrepant from our measurements of Na I and 
K I, respectively. This i s very consistent wi th the fact that HD 209458b seems to have 



a temperature inversion (IKnutson et al.l 120081 ). implying the need for an absorber such as 



TiO, VO, or something else — though as noted above, there is significant debate about which 
absorber (s) can properly account for inversions. 



^The targets for this survey were selected in 2007 (earlier in the case of HD 189733b) when far fewer 
transiting exoplanetary targets were known, and all four stars remain among the brightest known to have a 
transiting planet. 
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As we noted in ^ we are just now entering a time when we can do this sort of compar- 
ative atmospheric study for exoplanets. In our small sample here, we have noted two cases 
that superficially favor different models (HD 149026b and HD 209458). This drives home 
the importance of obtaining more data and creating data-driven models. 



5.6. Future Work 



The spectra that we have obtained are a very rich resource for additional studies. We 
are currently searching for Ha absorption in these four targets (Jensen et al. 2011b, in prepa- 
ration). As noted above, both HD 189733b and HD 209458b have had detections of Lya, ap- 
parently from the unbound p o rtions of their atmosphere, or exosphere. Using low-resolution 
ACS spectra, iBallester et al.l (120071 ) also detected absorption in HD 209458b corresponding 
to the B aimer edge of H This implicitly predicts that Ha absorpti on should occur in HD 
209458b's atmosphere, though past searches have turned up empty (jWinn et al.ll2004j ). A 
detection of Ha (n = 2 3) in any exoplanetary atmosphere would be a first. Combined 
with an existing Lya detection it would be a significant result, allowing for the calculation of 
a hydrogen excitation temperature and implications for the atmospheric physics. The many 
spectra that we have obtained over many epochs also present an interesting opportunity for 
work on stellar variability, including exploring activity-related intera ctions between th e stars 
and planets in these systems (e.g., the preliminary work reported by lArai et al.ll2010[ ). Most 
importantly, additional targets remain a possibility for this program. Observations of one 
such target, a hot Neptune-class planet, have been completed and are awaiting analysis. 



SUMMARY 



We have presented observations and analysis of the atmospheric transmission spectra 
of four transiting exoplanets using data taken with the HRS instrument on the HET. We 
searc h for Na I and K I ab sorption in four targets. Our target list includes HD 189733b, in 
which iRedfield et al.l (120081) detected Na I ab sorption. In an independent analysis of the same 



da taset, we confirm th e 



Redfield et al.l (120081) results. W e also tentatively confirm the results 



of ICharbonneau et all (I2OO2I 1 and ISnellen et al.l (120081 1 for Na I absorption in HD 209458b. 



Furthermore, we find a hint of Na I absorption in HD 149026b but cannot make a conclusive 



■^By the Balmer edge of H I we mean the combination of the n = 2 — > n = 00 transition edge at 3646 A, the 
associated bound- free absorption from the n — 2 state at nearby but shorter wavelengths, and bound-bound 
absorption from n = 2— >n3>2at nearby but longer wavelengths. 
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detection because of the large random and systematic errors. However, we suggest this is an 
excellent target for follow-up observations. We do not find evidence of Na I absorption in 
HD 147506b or K I absorption in any of our targets. Future work includes the analysis of a 
fifth target for which observations are complete, and searching for exospheric absorption in 
our targets. 
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Fig. 1. — Example of the spectra in the echelle orders containing Na I A5889, 5895 Ca I 
A6122 K I A7698 for HD 189733, along with the spectra of the telluric standard star HD 
196867. All spectra are reduced from the raw echelle images but are not processed with the 
additional steps we discuss in §2.21 (including removal of stellar and ISM features from the 
telluric spectra, telluric subtraction, and normalization). The target spectra are shown on 
the count scale of the y-axis, but the telluric spectra ars scaled to have the same maximum as 
the corresponding target spectra, and then arbitrarily shifted by 25% of that same maximum. 
The global shapes of the spectra are dominated by the blaze function of the echelle. In the 
central panel, note the lack of telluric features in the telluric standard spectrum near Ca I. 
Also note that the line of interest is closer to the peak of the blaze function than either 
Na I or K I. In the bottom panel, note the broad shallow feature near 7773 A. While this 
feature does not affect our measurements at K I, similar (but less dramatric) stellar and 
interstellar features are modeled and removed from the telluric spectrum near Na I and K I 
before subtraction from the primary science target. 
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Fig. 2. — The various transmission spectra at Na I A5889, 5895. Clockwise from upper 
left:— HD 147506, HD 149026, HD 209458, HD 189733. The binned transmission spectrum 
is shown in red, with the error bars of the unbinned spectrum shown in gray beneath the 
binned spectrum. For reference, the master out-of-transit spectrum of the star is shown at 
the top of each panel in green, scaled down to fit the plot, with a dotted green line showing 
the zero level and the top of each plot being unity in the normalized spectrum. A blue 
dashed line shows the zero point for the transmission spectrum. The purple vertical dot- 
dashed lines define the bandpass that is integrated to make our absorption measurements; 
from these lines to the edges of the plot define the comparison bandpasses. The right-hand 
y-axis shows the absorption in terms of relative apparent planetary radii. Values less than 
1.0 are not shown on this axis because they are unphysical (and should only occur due to 
noise or artifacts). 
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Fig. 3.— Same as Fig. 1 but for Ca I A6122. 
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Fig. 4.— Same as Fig. 1 but for K I A7698. 
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Fig. 5.— EMC analysis results for Na I A5889, 5895. Clockwise from upper left:— HD 147506, 
HD 149026, HD 209458, HD 189733. The distributions for the three versions of the EMC 
(see §3.2p are plotted in different colors, shown in the key: green is the "in-in" distribution, 
red is the "out-out" distribution, and blue is the "in-out" distribution. 5000 iterations are 
used in each case. The master "in-out" measurement is shown as a vertical dot-dashed line. 
Note that the scale of the x-axis is chosen to be consistent between this plot and the similar 
plots in Fig. [6] and Fig. [71 because the width of the distribution is the important quantity. 
Conversely, the scale of the y-axis is varied in each panel to show maximum detail. 
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Fig. 6.— Same as Fig. but for Ca I A6122. 
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Fig. 7.— Same as Fig. but for K I A7698. 
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Fig. 8. — A plot of all absorption results as a function of scale height of the atmosphere, 
with different color-coded symbols for each element. The scale height is calculated using 
Jupiter's mean molecular weight of 2.22 g/mol and the radiative equilibrium temperature of 
each planet. Na I is plotted at the proper scale height for each planet, while Ca I and K I 
are slightly offset for clarity, la error bars are shown. A dot-dashed line is also shown at 
(no absorption) for reference. 
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Fig. 9. — Same as Fig. |8] except with the transmission signal measurements converted into 
apparent planetary radius, relative to each planet. Note that the A in Ra is not a specific 
wavelength in this case but an integrated value based on the bandpasses defined in Table 
[31 A key to the symbols is shown on the plot; solid points are measured values, while open 
points are the model values discussed in the text, color-coded and aligned on the x-axis with 
the measurement it is intend to be compared to. 



Table 1. Observations 



Target 


In Transit'' 


Out of Transit* 


Ingress /Egress'' 


Telluric 


Observations 




# of Obs. 


Time (s)'' 


# of Obs. 


Time (s)'' 


# of Obs. 


Time (s)*" 




# of Obs. 


<Time> (s)*" 


HD 147506b"= 


28 


16508 


84 


49687 


5 


3000 


HD 156729^^ 


34 


567 


HD 149026b 


22 (1) 


12718 (600) 


90 


52734 


5 


3000 


HD 156729<^ 


34 


567 


HD 189733b 


38 (3) 


20140 (1610) 


173 (1) 


79011 (~0) 


15 (1) 


8213 (600) 


HD 196867 


45 


291 


HD 20<)ir,8l) 


29 {(>) 


17310 (3r)10) 


72 (3) 


1317.^) (1800) 


i 


2100 


IID 218015 


21 


80 



''Numbers shown arc all observations taken, while the numbers in parentheses indicate observations that are discarded for various reasons, 
primary cloud-induced low S/N, but also various reduction artifacts. The "In Transit" and ""Out of Transit" columns are mutually exclusive, 
defined by observation midpoint, while the "Ingress/Egress" observations are a subset of the "In Transit" observations, again defined by 
observation midpoint. For our Na I analysis, one additional 600 s out-of-transit observation is discarded for HD 189733b. For our Ca I analysis, 
one additional 443 s in-transit observation is discarded for HD 149026b. Finally, for HD 147506b, several observations are discarded for K L 3 
in-transit observations totaling 1800 s, and 6 out-of-transit observations totaling 3600 s. 

''Times for primary science target exposures are total times, while telluric exposure times are average times, as the telluric observations are 

always used independently and are never coadded. 

<:HD 147506b is also known as HAT-P-2. 

"^HD 147506b and HD 149026b use the same telluric. The values in the table this telluric represent the same observations, i.e., the 34 
observations listed do not distinguish between which primary science target was originally associated with the telluric observation. 



Table 2. System and Orbital Parameters of Targets 



Target 


Reference Epoch 
(HJD-2400000) 


Period 

(days) 


Transit Duration 
(minutes) 


a/R, 


Rp/R, 


Inclination 

(°) 


References 


HD 147506b 
HD 149026b 
HD 189733b 
HD 209458b 


54387.4937 ± 0.00074 
54456.7876 ± 0.00014 
54279.43671 ±0.000015 
52854.825 ± 0.000135 


5.633473 ± 0.0000061 
2.87589 ±0.00001 
2.2185757 ± 0.00000015 
3.5247455 ± 0.00000018 


257.76 ± 1.87 
194.4 ±8.7 

106.1 ±0.7 

184.2 ± 1.6 


8.99lOi? 

7 11+0.03 

8.81 ± 0.06 
8.76 ± 0.04 


0.07227 ± 0.00061 

n 0401+0 0018 
0.15436 ±0.00022 
0.12086 ± 0.0001 


86.7t;i? 

85.58 ±0.06 
86.71 ±0.05 


1, 2, 3, 4 
3, 5, 6 

3, 7, 8, 9 
3, 10 



Ref erenc es. — (1*) |PaI et al.l boidb. (2) iBakos et all boidb. (3')lTorres et aO ll2008h. (4) IWinii et a.1.] ||2007| '). (5) ICarter et al] | |2009| '|. (6) IWolf et al.l 
j2007l ). (7) 'Agol ct al.' ( ^OldV (Sl lTriaud et al.l l l200g^ . Ol lBouchv et al.l | |2005|) . (lO)'Winn ct al.' ('2005'). This table also made use of the Exoplanct Data 
Explorer ( http:/ /exoplanets.orgj ) and the Exoplanet Transit Database (http://var2.astro.cz/ETD/indcx.php). 
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Table 3. 


Absorption Line Data 


Element 


Wavelength(s) (A) 


Bandwidth (A)^ Use 


Nal 


5889.9510, 5895.9242 


12 Alkali metal resonance line 


K I 


7698.9645 


8 Alkali metal resonance line 


Cal 


6122.2217 


8 Control line 



'^The bandwidth over which integration is performed. The band is centered at 
the wavelength in the first column, or in the case of multiple lines for Na I, their 
average. We also use reference bands of equal bandwidth immediately to the blue 
and red of these bands. 



Table 4. Absorption Results 



Star Nal Cal KI 

HD 147506b (1.88 ± 1.18) x 10"^ (-0.078 ± 0.74) x 10"'' (2.17 ± 1.02) x 10"'' 

HD 149026b (-3.16 ± 2.06) x 10'"^ (1.15 ± 0.78) x 10"'' (-0.71 ± 1.09) x 10"^ 

HD 189733b (-5.26 ± 1.69) x 10"^ (0.72 ± 0.63) x 10"^ (0.77 ± 1.04) x 10"^ 

HD 209458b (-2.63 ± 0.62) x lO"'' (-0.36 ± 0.36) x 10"^ (-0.89 ± 0.67) x 10"^ 



